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Abstract 

We examine prospects for detecting the neutral Higgs bosons of minimal 
supersymmetric models (MSSM) when their decays into neutralino pairs are 
kinematically allowed. The best signature appears to be Hh,Hp — > Z2Z2 
41 + I^j,. We argue that Standard Model contributions to this signature are 
negligible, and examine regions of MSSM parameter space where the four 
lepton mode should be observable at the Large Hadron Collider. The same 
signal can also come from continuum neutralino pair production. We pro- 
pose a set of cuts to illustrate that the neutralino decay mode of the Higgs 
bosons provides a viable signal over a substantial range of model parameters, 
and show that it may be separable from continuum neutralino production if 
sufficient integrated luminosity can be accumulated. 
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Typeset using REVT^ 



I. INTRODUCTION 



The derivation of large radiative corrections to Higgs boson masses and couplings in 
the Minimal Supersymmetric Standard Model (MSSM) has led a number of groups 
to re-evaluate prospects for SUSY Higgs detection at various colliding beam facilities. Most 
of these studies have focused on regions of MSSM parameter space where various standard 
model (SM) decay modes {e.g. 77 and ZZ or ZZ* — >• 4£ ) of the SUSY Higgs bosons are 
detectable above background; parameter space choices were selected such that SUSY particle 
masses were large so that Higgs boson decays to sparticles were kinematically forbidden. A 
region of parameter space roughly spanning pseudoscalar Higgs mass rriHp ~ 100 — 300 GeV 
and ratio of Higgs vevs tan (3 ~ 4 — 10 (for tan/9 > 10, the observability of the signal is 
sensitive to the detectability of the rf decay modes of the Higgs bosons [0) was found where 
none of the SM decay modes were detectable [|10|- Very recently, it has been argued [|ll|] that 
Higgs bosons, produced in association with t-quark pairs and identified via their dominant 
bb decays may fill this "hole" , provided that sufficient 6-tagging capability can be achieved. 

Over the last year or two, several groups 0] have studied grand unified models within 
the supergravity framework and have shown that it is quite possible to construct models 
consistent with constraints from colliders, proton decay experiments and cosmology. Inter- 
estingly, these analyses generally find that the sparticles are all considerably lighter than 1 
TeV, and further, that the lighter chargino (Wi) and the two lighter neutralinos {Zi and 
Z2) frequently have masses in the range 50-150 GeV so that these may well be accessible 
via the decays of the MSSM Higgs bosons. 

In Ref. 0, we studied how these supersymmetric decay modes would affect the phe- 
nomenology of the Higgs sector. We showed that the chargino and neutralino decays of 
MSSM Higgs bosons have substantial branching fractions when kinematically allowed, and 
can sometimes even dominate the usual decay modes. Decays to top squark pairs can sim- 
ilarly be significant since ti, the lighter of the two t-squarks, can be much lighter than all 
other squarks. These new decay channels lead to a diminution of the rate into the stan- 
dard decay modes so that the above-mentioned "hole" in parameter space becomes larger. 
Moreover, the region of parameter space where more than one of the Higgs bosons is visi- 
ble above background (leading to unambiguous evidence for a non-minimal Higgs sector) is 
substantially diminished. 

One may also ask if the new supersymmetric decay modes of Higgs bosons can lead to 
new avenues for detection. Decay modes of the neutral Higgs bosons into titi and into 
WiWi will lead in general to final states containing — 2 leptons plus jets plus I^j.; at 
the proposed CERN Large Hadron Collider (LHC), such signals are expected to be buried 
beneath standard model backgrounds from processes such as top quark pair production. A 
more promising signature may be found by searching for Hp,Hh Z2Z2, where the 
neutralinos decay leptonically via Z2 iiZi. Such a process leads to a final state with as 
many as four isolated leptons plus I^j., which is expected to have small SM backgrounds. In 
addition, kinematic information from the 4i + final state can yield information on not 
only the Higgs boson masses, but also on the masses of Z2 and Zi |^. Signals for invisible 
decays of Higgs bosons via He Z\Z\ have also been considered p. 

In this paper, we seek to expand upon and improve the calculations presented in Ref. |0] 
regarding signals from supersymmetric decays of the Higgs bosons. In particular. 
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we incorporate radiative corrections to the SUSY Higgs boson masses and couplings, 
from both top and bottom Yukawa interactions, and effects from mixing between third 
generation squarks, 

we include Higgs boson decays into squarks and sleptons — these can be very important 
since decays to ti pairs may be allowed, 

we have included Higgs production via bb fusion, which is the main production mech- 
anism for large values of tan (3, 

we have included effects from non-degenerate squark and slepton masses, which can 
lead to enhanced |jl3| leptonic decays of the neutralino, resulting in a much larger 



signal. 

The rest of this paper is organized as follows. In Sec. 2, we describe in some detail 
the improvements and extensions we have made to our earlier calculation. In Sec. 3, we 
present cross-sections for the various signals as a function of model parameters. In Sec. 
4, we present background calculations, and suggest a set of experimental cuts useful for 
extracting the signal. We also estimate efficiencies due to these cuts. In Sec. 5, we present 
our conclusions regarding detectability of the SUSY Higgs bosons via both the standard 
model decays as well as via the supersymmetric Ai signal, along with a summary of our 
results. 



II. CALCULATIONAL DETAILS 

As in Ref. [0], we have computed the masses and mixing angles in the Higgs boson sector 
using the 1-loop effective potential. We have, however, improved on our previous calculation 
in several respects. We now include the effects of both top and bottom family Yukawa 
interactions (the latter interactions, which were neglected in Ref. 0, can be important if 
tan /3 is large) as well as effects from supersymmetric and SUSY-breaking trilinear scalar 
interactions of the scalar top (scalar bottom) Higgs system. These trilinear couplings lead 
to mixing between the L- and R- type sfermions, and so further reduce the mass of the 
lighter mass eigenstate. This can be especially important for t-squarks, whose mass is 
already expected to be smaller than that of other squarks on account of the large top 
Yukawa coupling, since the branching fraction for the titi decay mode can be very large if 
< rrit. Finally, we have also included D-term contributions to the effective potential, 
but generally speaking, their effects are small. 

These modifications lead to improvements in the calculation of some of the decay modes 
of the Higgs bosons. First, radiative corrections at the same level as those mentioned in the 
last paragraph are incorporated into the calculation of the Hh-Hi-Hi vertex. Second, Higgs 
decays into sfermions are now added. As mentioned, the sfermion masses including D-terms 
have been calculated with intra-fiavor mixing from non-zero values for At, A^, and fi. These 
modifications to squark masses and mixings also affect the loop decays of the Higgs bosons 
into two gluons which, in turn, affects Higgs production by gluon fusion. 

In the MSSM, the bottom quark Yukawa coupling is inversely proportional to the pa- 
rameter cos (3. Hence, the H bb width (here, if is a generic Higgs boson) as well as the 
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subprocess cross section a{bb H) are enhanced in regions of large tan/5. In our calcula- 
tions, we include Higgs boson production cross sections via gg and bb fusion. The formula 
for a{gg — H) is well known |T5[; for bb fusion, we use the result 



a{pp ^bb^ HX) 



167r2 1 1 



T{H bb) 



-tFdw X 



4:9 Xh[l^rnl/mjj,ml/mjj) 
[' -( Q^) A/,(x/r, Q') + (6 ^ b)) 

J T X 



where r = mjj/s, Dh/p{x, Q'^) is the b parton distribution function in the proton, and 
Fdw = [16.2 - 4.28 ln(m^^) + 0.31(ln(mjf))2]-i 



(2.1) 
(2.2) 

(2.3) 



(with tuh in GeV) is a fit to the results of Ref. 0] to incorporate contributions of higher 
order graphs to the bb fusion mechanism. We use the EHLQ Set 1 [O] b and g parton 
distributions, and take the SM parameters to be = 5 GeV, sin^ = 0.23, ^^(M^^) = 
0.118 with A^iQCD) = 0.177. 



It has been pointed out in Ref. [|I^ that leptonic decays of neutralinos can be enhanced by 
large factors if squarks are significantly heavier than sleptons (as is the case in the "no-scale" 

limit), and the ZZ^Z^ coupling is dynamically suppressed. This suppression is common 
in supergravity models with radiative electroweak symmetry breaking since, in this case, 

~ rag P,P^ — as a result, Z\ and Z^ are respectively mainly U(l) and SU(2) gauginos 
and so have suppressed couplings to the Z boson. We calculate slepton masses as usual by 
using renormalization group equations to evolve sfermion masses from a common GUT scale 
scalar mass to their weak scale values |1^ . In our computation of the neutralino mass matrix 



we have, as usual, assumed that the MS gaugino masses unify at some ultrahigh energy 
scale. We then evolve these down to the electroweak scale. However, in the following, 
we present our results in terms of the physical pole gluino mass which we relate to the 

Finally, we note that we have not 



corresponding MS mass using the result of Ref. 
included QCD corrections to Higgs production and decay via loops [p^ 



III. THE 4-LEPTON SIGNAL FROM NEUTRALINO DECAYS OF HIGGS 

BOSONS 

Once the various Higgs boson production cross sections via gg and 66 fusion as well 
as the H Z^Z^ and Z2 — * llZx branching ratios are known, the total rate for cr(pp — »• 
H — »• Z2Z2 — >■ 4£ + 2Zi) can be calculated. The cross section is a function of the MSSM 
model parameters (m^, m^, m^^, m^^, rrij^^, /i, tan /3, m//p, A^, A^, m^); the slepton masses are 
related as given in Ref. [|T^]. Of course, in supergravity models with radiative electroweak 
symmetry breaking, there is some correlation amongst these parameters: typically, /i scales 
with m^, and vhh^ is strongly correlated with /i and the universal GUT scale scalar mass. 
The parameter choices used in this paper to illustrate our results are inspired, but not ruled, 
by the supergravity mass relations. 

We begin by showing in Fig. 1 a contour plot in the vnn^ vs. tan/5 plane of aipp H ^ 
Z2Z2 4£ + (£ = e or /x) in /6, at = 14 TeV. We take = -/x = = 300 GeV, 
while At = Ah = 0, and rrit = 165 GeV. The scalar top masses are set at their default values: 
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m? = m? — 50 GeV^, and m? = m? — 100 GeV^. The reEfion in black is excluded by the 
non-observation of supersymmetry signals in experiments at LEP as parametrized in Ref. 
but taking into account the recent limit of 63.5 GeV on the mass of the SM Higgs boson 
P^ . Fig. la shows results for H = Hp, while Fig. lb shows results for H = H^. We note 
the following features. 

• The Hp — >■ M cross section exceeds 500 fh for small values of tan f3 and mn^, ~ 300 GeV, 
just below threshold for Hp — » ti decay whereas the corresponding cross section from 
Hh does not exceed 100 fh. Over a wide range of parameters, away from tan/? = 1, 
the two processes give comparable cross sections. The large difference between the 
scalar and pseudoscalar contributions to the signal for small values of tan/3 mainly 
comes from the difference in their SUSY branching fractions 0]. 

• The signal from the decays of Hp and H^ are separately larger than 5 fh for 200 GeV 
< rriH^ < 400 GeV. Since H^ and Hp are expected to be roughly degenerate over the 
range of masses where the signal is significant, this corresponds to a total of 250-25K 
4i events at the LHC before any selection cuts, assuming a data sample with 50 fh^^. 
Interestingly, the signal has a larger rate in the region below the dashed line, where 
the lighter chargino is heavier than 90 GeV, which we take to roughly represent the 
supersymmetry reach of LEP II (We note that this region is rather sensitive to the 
precise value of the chargino mass reach that is assumed). The Higgs bosons Hi and 
Hp may themselves be directly accessible at LEP II. The range of parameters where 
this is possible is discussed in the last section and illustrated in Fig. 7 for one choice 
of parameters. 

• Although we have not shown this, we have checked that for values of tan [3^4, both 
Hh and Hp cross sections are dominated by gluon fusion, while hh fusion, which was 
ignored in Ref. |^, dominates for tan/? ^ 10. 

• The cross section for M + I^j, events has an observable rate in part of the "hole region" , 
where 3 < tan j3 < 10—20 and 100 ^ ~ 200—300 GeV. This is the region where no 
SM decays of MSSM Higgs are observable (except for possible observation of H —>■ hh). 
Of course, the range over which this signal might be observable is sensitive to fi and rrig 
since the cross section drops sharply to zero when the kinematic limit for the decays 
Hp h Z2Z2 is approached. 

The contours of the 4i cross section are shown in the fi — tan (3 plane in Fig. 2 for (a) Hp 
decays and (6) H^ decays. The pseudoscalar mass is fixed to be 250 GeV, which is within the 
hole region for the light Wi and Zi^2 case illustrated here. Other parameters are as in Fig. 
1. We see that the signal cross section exceeds 10 fh for a wide region of parameter space. 
Again, the black region is excluded by LEP constraints and the dashed line is the contour 
= 90 GeV (the region above the dashed contour has < 90 GeV). In the upper 
corners of Fig. 2, rrif < 45 GeV, because of r Yukawa interactions which have substantial 
effects for large values of tan/? |^3[. While the signal is small for small values of \fi\, it is 



instructive to see that the 4£ rate appears to be observable for \fi\ ~ rUg as expected in 
supergravity models with radiative EW symmetry breaking. Finally, we observe that the 
pseudoscalar Higgs boson tends to give a larger signal than the heavy scalar, and further, 
that the signal is significantly larger if /i happens to be negative in our convention. 
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In Fig. 3, we study the dependence of the signal on the squark mass for a fixed value 
of rrig. Here, is fixed at 250 GeV; other parameters are as in Fig. 1. We see that 
the 4i cross section rapidly drops off as rrig increases from rrig to larger values. For values 
of tan/3 ^ 5, the cross section shows a slow increase for rather large squark masses. As 
expected, the squark mass dependence of the cross section is mostly determined by that 
of the leptonic branching fraction of the Z2. For rriq ~ rrig, the sleptons are considerably 
lighter than the squarks resulting [O] in an enhancement of Z2 iiZi decays, and hence, 
in the 4i cross section. The increase in the cross section for large values of rrig comes from 
interference between various amplitudes. We see that the signal has an observable rate even 
if squarks are significantly heavier than gluinos. Amusingly, the cross section is largest below 
the dashed line where the chargino is likely to be undetectable at LEP II. 

Before closing this discussion we mention that we have checked that the 4i cross section 
is relatively insensitive to the A-parameters or to the precise value of except when the 
threshold for a new decay is crossed. This is especially important when considering the 
variation of the signal with At, since H titi decays become kinematically accessible when 
\At\ becomes very large. Finally, we have also studied the variation of the signal with rUg. 
For rriHp = 250 GeV, and —fi = rrig = rrig, the signal exceeds 5 fb for gluinos as heavy as 
350-400 GeV. This is reasonable since then, m^^ ~= rrig/ 3 is approaching the kinematic 

boundary of Hp Z2Z2 decays. 



IV. BACKGROUNDS, CUTS AND EFFICIENCIES 



We have seen that the neutralino decays of the pseudoscalar and the heavy neutral scalar 
Higgs bosons of the MSSM can result in several hundred to several thousand events with 
four isolated leptons together in an LHC data sample of 50 fb~^. These events should be 
very distinctive since, except for QCD radiation, they would be free of hadronic activity. 
Within the SM, isolated leptons mainly arise from the the production and decays of the W 
and Z bosons and of the top quark. The main background from ZZ pair production can be 
efficiently removed by vetoing events where like-flavour, opposite sign leptons reconstruct the 
Z mass within ±10 GeV. We have verifled this by generating lOK ZZ events using ISAJET 
(and forcing the leptonic decay of the Z bosons). No events are found after the mass cut, 
resulting in an upper bound on the 4£ cross section of about 0.003 fb. This would, of course, 
also remove the signal if the decay Z2 Z^Z is kinematically accessible; fortunately for the 
range of parameters where we flnd the signal to be substantial, this decay is kinematically 
inaccessible, so that this requirement results in very little loss of the signal. The signal can 
also be mimicked by ttZ or 4t production. These backgrounds can easily be removed by 
vetoing events with a central jet in addition to the Z-mass veto already mentioned. The 
main SM physics background thus comes from electroweak multi-VT production. At the 
LHC, the trilepton cross section from ?)W production has been shown |2^ to be about 2 



fb, so that the background to the signal from 41^ production should be negligible. Because 
the cross section for production at the LHC is very large, one may also worry that these 
may provide a signiflcant background when the leptons from the daughter bottom quarks 
are accidently isolated; a simulation of 300K ti events with forced top quark decays yielded 
no background, giving a limit cr(tt) < 4 fb. We have been unable to identify any signiflcant 
SM sources of physics backgrounds to the SUSY Higgs boson signal. 
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Within the MSSM framework, however, the continuum production of pairs ||2^ can 
also resuh in the same signal. These are produced by qq annihilation via s-channel Z 
exchange or t- and w-channel squark exchange. While the detection of these continuum 
neutralino pairs would in itself be very exciting, it is interesting to ask whether Z2Z2 pro- 
duction via Hh and Hp decays is distinguishable from the continuum production of Z2 which, 
in effect, is the background to our Higgs signal. Toward this end, we have shown in Fig. 4 
the cross section for 4£ production via continuum Z2Z2 production as a function of tan/?. 
To compare with Z2 production via Higgs boson decays discussed in the last section, we 
have taken m-g = ±/i = 300 GeV and illustrated our results for nig = nig (solid line) and 
niq = 2nig (dashed line). The pseudoscalar Higgs boson mass is fixed to be 250 GeV and 
other parameters are fixed as in Fig. 1-3. We see that while the cross section is sensitive to 
the squark mass it is relatively insensitive to tan/3 over the region where the cross section 
is significant. Furthermore, for nig = rriq, the cross section is 10-30 fb which is generally 
comparable to, or smaller than, the cross sections in Fig. 2 for a wide range of parameters. 

We use IS A JET 7.07 to simulate the M Higgs boson signal. Since explicit SUSY 



Higgs production has not yet been incorporated into this code, we simulate the production 
of Hh or Hp by decaying the SM Higgs scalar into a Z2 pair and forcing the SUSY decay 
mode; the total cross section is then normalized to the results of Fig. 1-3. 

We use the toy calorimeter simulation package ISAPLT to model detector effects. We 
simulate calorimetry with cell size At] x A(j) = 0.1 x 0.1, which extends between —5 < 1] < 5 
in pseudorapidity. We take electromagnetic energy resolution to be 10%/a/Et © 0.01, 
while hadronic resolution is © 0.03 for |?7| < 3, and 100%/v^ © 0.07 for 

3 < \ri\ < 5, where © denotes addition in quadrature. Jets are coalesced within cones 
of i? = y/ Arj'^ + A(j)'^ = 0.7 using the ISAJET routine GETJET. Hadronic clusters with 
Et > 50 GeV are labelled as jets. Muons and electrons are classified as isolated if they 
have pt > 10 GeV, \'ri{i)\ < 2.5, and the visible activity within a cone of i? = 0.3 about the 
lepton direction is less than Exlcone) = 2 GeV. 

We then impose the following cuts designed to select signal events, while vetoing SM 
backgrounds from ZZ and tt production: 

• require two isolated leptons with pt(^) > 20 GeV to trigger the event, 

• require two more isolated leptons with pt{^) > 10 GeV, 

• require all opposite sign but same flavor dilepton pairs to have invariant mass ni{ii) < 
80 GeV or ni{U) > 100 GeV, 

• require number of jets n(jets) = 0. 

A cut on I^j, could be considered instead of the above dilepton mass cut. However, the 
J^rp spectrum from the signal is not so hard, while forward jet production and energy mis- 
measurement in ZZ events can lead to substantial so that the dilepton mass cut wins 
over an cut in rejecting background while preserving signal. 

We should stress that there are no significant SM backgrounds even before the last cut. 
Backgrounds involving the Z are efficiently removed by the ni{ii) cut. The main physics 
background would then be expected to come from 4t production which gives a 4i cross 
section around 0.05/6 at the LHC, even before acceptance cuts The last cut which 
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removes about 40% (70%) of the signal for a Higgs mass around 200 GeV (400 GeV) has 
been imposed to separate neutrahno production from the cascade decays of gluinos and 
squarks which can also produce multilepton events at observable rates. 

In order to give the reader an idea of the impact of the cuts on the 4i signal, we have 
shown these cross sections for illustrative choices of rriHp and tan j3 in Table I. We have added 
the contributions from Hh and Hp decays since these are expected to lead to kinematically 
similar events. We have also shown the continuum background for the same choices of 
parameters. The following points are worth noting. 

• For the choice of parameters in the Table, the signal exceeds the background for 

up to somewhat beyond 300 GeV, where ti decays become accessible. Also, the cross 
section corresponds to an event rate 100-1000 events after cuts in a 50 fb~^ data 
sample, compared to a continuum background of 50-100 events. 

• The signal efficiency varies between 5 and 10 percent depending on the model param- 
eters; Fig. 1-3 can thus be used to estimate the signal after the cuts. 

It should, of course, be kept in mind that the SUSY parameters are not known, and 
that the total background (and signal) rate could be considerably different (even for similar 
values of tu-^^) from our estimate in the Table. An excess of the 4i events relative to the 
background in Table I would, therefore, not enable us to cleanly infer a Higgs boson signal. 
Instead, we consider the possibihty of separating the signal from the continuum background 
by studying the invariant mass distribution of the four leptons: for the signal, we must have 
m(4£)< ruH — 2m^^, while the background should exhibit a rather broad distribution. Since 
these distributions are determined by the Higgs boson and neutralino masses, we expect 
them to be relatively insensitive to variations in model parameters which result in similar 
values of mn, and w^^^- 

Toward this end, we have shown in Fig. 5 and Fig. 6 these distributions for the signal 
plus background (solid) and the Z2Z2 continuum "background" (dashed) for the six cases 
in Table I. For the smaller values of mupi the solid histograms are dominated by the signal 
and differ considerably from the dashed background histograms (note that these are shown 
using a log scale). As anticipated above, the solid and dashed lines indeed coincide for 
m(4£) > niHp — 2m^^. In order to decide whether the solid and dashed histograms are 
indeed distinguishable, we have computed the total fo^' the difference between these two 
histograms, normalized to the same number of events, with the event number given by the 
signal plus background cross section times an integrated luminosity of 50 ph~^ . In this 
computation, we have used twelve 20 GeV bins between m(4£) = 60 GeV and m(4£) = 300 
GeV. In our simulation, we have about 1100, 900 and 600 events for = 200, 300 and 
400 GeV (after cuts) for the tan P = 2 case, and about 550 events for each value of niHp for 
tan j3 = 10. The resulting total is shown in the last column of Table I. It thus appears that 
for mfjp ^ 300 GeV, the distributions are sufficiently different that the solid line is unlikely 
to be a chance fluctuation of the continuum background (for > 26.2, this probability is 
smaller than 1%). Some remarks are, however, in order: 

• Despite the fact that we have normalized the solid and histograms to have the same 
number of events, our conclusion clearly depends on the relative number of Higgs 
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initiated and continuum Z2Z2 events. Here, we have used the rate as given by the 
MSSM for the values of parameters motivated by supergravity models. 

• For the first two cases in Table I, the number of events in our simulation is comparable 
to the expected number in a data sample of 50 fh~^. We thus expect our computation 
of the total to be a reasonable reflection of the experimental situation. For the 
tan/5=10 case, with rau^ = 200 GeV (300 GeV), the signal cross section is much 
smaller so that about 200 fb~^ (400 fb~^) of integrated luminosity need to be collected 
in order that the fluctuations in our simulation are of comparable magnitude to those 
in the data. We thus conclude that while 50 fb~^ of data may suffice to enable one 
to distinguish the Higgs signal from continuum Z2Z2 production for smaller values of 
tan j3, integrated luminosities of 200-400 fb~^ may be necessary if tan (3 is large. 

• We should bear in mind that we have used only the shape of the 4i mass distribution to 
try to untangle the Higgs signal from continuum Z2Z2 production without any regard 
for rate or other event shape variables. It would be interesting to explore whether 
other distributions serve as better discriminators of the Higgs from the continuum 
background. 



V. CONCLUDING REMARKS 

If low energy supersymmetry is to provide a rationale for the stability of the electroweak 
scale, sparticles must all be lighter than about 1 TeV. In models where gaugino masses are 
unified at an ultra-high unification scale, the lighter chargino and the two lightest neutralinos 
are frequently lighter than 100-150 GeV, and are thus expected to be kinematically accessible 
in the decays of the heavier Higgs bosons of the model. In a previous paper |]^, we had shown 
within the MSSM framework that these SUSY modes dominated the decays of the heavier 
neutral Higgs boson, and particularly, the pseudoscalar Higgs boson over a wide range of 
SUSY parameters. This has two important consequences. The down side is that it reduces 
the cross sections for the 77 and ZZ or ZZ* modes which form the usual signal for Higgs 
bosons at hadron colliders. These SUSY decays, however, open up new possibilities for Higgs 
boson detection, the most promising of which is the 4i signal from the Z2Z2 decays of Hp or 
Hh- Here, we have improved on our previous computation of this signal, and also explored 
how it varies with model parameters. In this connection, we have used supergravity models 
as a guide in order to restrict the parameter space, and make this exploration tractable. 

For the convenience of the reader, and also because the projected energy of the LHC has 
been reduced to 14 TeV, we first briefiy review the detectability of the SM decay modes of 
the Higgs boson at the LHC. For the 77 signal, we require the center of mass scattering angle 



satisfy cos 6** < 0.8 [^]. As before 0, the background is assumed to come from continuum 
pair production via qq,gg 77 where the photons have m^^ = rriH ± 1%. For the ZZ 
and ZZ* signal, we have required that the four leptons reconstruct to the Higgs boson mass 
within a mass- dependent resolution given in Ref. We have considered the ZZ* — > 4i 
signal only for mn > 130 GeV, in which case backgrounds from continuum ZZ* and Z7* 
production have been shown to be negligible ||30[| . 



The regions of the rriHp — tan /5 plane where these various signals are observable using 
the 99% CL criterion described in Ref. is illustrated in Fig. 7 for (a) rrig = rrig = —fi = 1 
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TeV, and (6) nig = rriq = —jj, = —300 GeV. In this figure, we have assumed At = Af, = 
and taken the integrated luminosity to be 50 fb~^. The legends in this figure appear on the 
side of the boundary where the signal is observable. The black region denotes the range of 



parameters excluded by experiments at LEP p2[ while the regions below the lines marked 



LEP190 and LEP175 can be probed in experiments at LEP with optimistic {a{H^Hp) or 
a{HeZ) > 0.05 pb, ^ = 190 GeV) and pessimistic ia{HiHp) or a{HiZ) > 0.2 pb, y/s = 175 
GeV) scenarios for performance of the LEP collider in its next phase. We note that: 

• We have used an integrated luminosity of 50 fb~^ because at the reduced energy, we 
found that the Hi — > 77 signal was essentially unobservable over the whole plane in 
Fig. 7b with a data sample of "just" 30 fb~^. This is a refiection of the well-known 
fact that the position of this contour is extremely sensitive to the assumptions about 
the detector. 

• The shaded region is where none of the neutral Higgs bosons of the MSSM are de- 
tectable at either the LHC or at LEP II, at least via the signals discussed above. As 
noted in the Introduction, it may be possible to detect Higgs boson signals even for 
parameters inside this "hole" if Higgs decays to r-leptons or bottom quark pairs 
PH are identifiable. 



• It is amusing to see that the shaded region actually shrinks in Fig. 7b. This is 
somewhat misleading because this shrinkage is due to the upward movement of the 
LEP190 curve. Conservative assumptions about the performance of LEP considerably 
increase the region where there is no signal either at LEP II, or at the LHC. We have 
also checked that increasing the value of the A-parameter lowers the LEP 190 curve. 
We have traced this to an increase in rriHi, corresponding suppression of the 
cross section for ZHe production: for At = 400 GeV or At = -700 GeV, the LEP 190 
curve roughly follows the LEP175 curve in the figure. Thus, the "hole" in Fig. 7b may 
well be considerably bigger than indicated by the shaded region even with optimistic 
assumptions about the performance of LEP II. 

• We note here that the main reason for the change in the 77 boundary when 
the parameters are altered from 1 TeV to 300 GeV is not the opening of the SUSY 
decays of Hi. The shift occurs primarily because the Higgs mass is altered (due to the 
incorporation of the improved radiative corrections) resulting in a different size of the 
background. 

• Finally, we note that as anticipated, the region of parameters where the signal from 
two different Higgs bosons is simultaneously detectable is greatly reduced for the case 
in Fig. 7b 

Turning to the neutralino signal for Higgs bosons, the cross sections for the 4i signal 
from Z2Z2 decays of the Higgs bosons are summarized in Figs. 1-3. We see that these cross 
sections can be as large as 500 fb for ranges of model parameters allowed by all known 
experimental data. It is also worth noting that, before experimental cuts, the signal exceeds 
10 fb over a large region of parameter space where there may be no visible SUSY signal 
even at LEP II. We have argued that there are no significant SM backgrounds to this signal. 
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Thus while the detection of four lepton events at such rates will be a signal for new physics, 
its identification as a Higgs boson signal requires that it be separable from continuum 
production for which the cross section is shown in Fig. 4, for /i = ^rrig as expected in 
supergravity models. We see that over a large range of parameters, this background is 
significantly smaller than the signal in Fig. 1-3. The efficiency with which this signal may 
be detected is shown in Table I, for cuts typical of LHC detectors. We see that this is 
typically 5-10%, so that the cross sections in Fig. 1-3 correspond to 25-1000 events in an 
experimental data sample of 50 fb~^. 

In order to assess whether the Higgs signal could be distinguished above the neutralino 
continuum, we examined the mass distribution of the four leptons produced via the decay 
of the Higgs boson. For illustrative values of model parameters, we showed that the shape 
of this distribution may serve to discriminate the Higgs signal from continuum neutralino 
production provided ^ 2mt. We have argued that for low values of tan P an integrated 
luminosity of 50 fb~^ is sufficient for this discrimination, but an integrated luminosity of 200- 
400 fb^^ is required if tan/5 = 10. We have, respectively, denoted these cases by crosses and 
open circles in Fig. 7b. As we can see SUSY decays of Higgs bosons are indeed detectable 
well into the "hole" region if LEP II is operated at about 175 GeV. Although this is not 
obvious from the figure, this may also be the case with optimistic assumptions about the 
performance of LEP II, since, as we mentioned the LEP II observability curve essentially 
follows the LEP 175 curve if \At\ is large, while the Higgs signal, is relatively insensitive to 
variations in At. 

In summary, we have shown that the processes Hh,p — >■ Z2Z2 M + I^j, lead to an 
observable rate for gold-plated four lepton events at the LHC for a significant range of 
SUSY model parameters. The parameter space region where the neutralino decays of SUSY 
Higgs bosons ought to be observable can be summarized as 

• 2mj > rriHp > 2^«^^ ~ 2mg/3, 

• niq ~ rrig so that << nig and ~ nig < 500 GeV 

• niHp ~ 200 - 350 GeV, 

We note that if SUSY parameters are in this region, there will also be a plethora of other 
signals via which SUSY will be detectable at the LHC p8| ). 

We have been unable to identify any significant SM background to this distinct SUSY 
Higgs boson signal. We have shown that this signal may be distinguished from continuum 
Z2Z2 production provided a sufficiently large integrated luminosity is available. Finally, this 
process may provide the only way to identify any Higgs boson of the MSSM if the model 
parameters are in the "hole" region (unless identification of their r and bottom quark decays 
turns out to be feasible), and further, that it may well provide the only identifiable signal 
for the notoriously hard to detect pseudoscalar Higgs boson. 
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TABLES 



TABLE I. Cross sections in fb at LHC for 4£ + ^j, events from supersymmetric processes. We 
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FIGURES 



FIG. 1. Contour plot of cross section in fb for a) pp ^ Hp ^ Z2Z2 M + ^j, and h) 
pp — Hh — Z2Z2 41 + ^j, events, in the mn^ vs. tan/3 plane, at \/s = 14 TeV. We take 
rrig = rUg = = 300 GeV, nit = 165 GeV, and At = Ab = 0. We also take = m| - 50 GcV^ 
and m? = m? — 100 GeV^. The region above the dashed lines corresponds to m~ < 90 GcV, the 
approximate reach of LEP II. The reach of LEP II for detection of Higgs boson signals is shown in 
Fig. 7b 

FIG. 2. Contour plot of cross section in fh for a) pp Hp Z2Z2 M + Ifrp and h) 
pp — Hh — Z2Z2 M + ^j, events, in the /x vs. tan (3 plane, at y/s = 14 TeV. Parameters are as 
in Fig. 1, except tuHj, = 250 GeV. The region above the dashed contour corresponds to < 90 
GeV. 

FIG. 3. Contour plot of cross section in fb for a) pp Hp Z2Z2 AC. + Ifj, and b) 
pp Hfi Z2Z2 AC + I^rj, events, in the rriq vs. tan fi plane, at = 14 TeV. Parameters are as 
in Fig. 1, except = 250 GeV. The region above the dashed contour corresponds to fn^^ < 90 
GeV. 

FIG. 4. Cross section in pb for continuum pp Z2Z2 M + ^j, production versus tan/3 for 
mg = ±iJ, = 300 GeV, at = 14 TeV. 

FIG. 5. DistribTition in rn{4:i) after cuts from signal plus background (solid), and background 
(dashes) for pp Hp, H^ Z2Z2 M + Ifj^ production for tan (3 = 2 for a) ruHp = 200 GeV, b) 
ruHp = 300 GeV, and c) ruHp = 400 GeV. Other parameters are as in Fig. 1. 

FIG. 6. Distribution in m(4^) after cuts from signal plus background (solid), and background 
(dashes) after cuts for pp Hp, Hh Z2Z2 M+^j, production for tan /? = 10 for a) niHp = 200 
GeV, b) ruHp = 300 GeV, and c) niHp = 400 GeV. Other parameters are as in Fig. 1. 

FIG. 7. Plot of discovery regions in the rUHp vs. tan/3 plane, at = 14 TeV, assuming an 
integrated luminosity of 50 fb~^. In a), we take rriq = nig = —jj, = 1000 GeV, while in 6), we take 
niq = nig = —jj. = 300 GeV, so that SUSY decay of Higgs bosons are allowed. Other parameters 
are as in Fig. 1. The Higgs boson reach of LEP II is sensitive to the value of the ^-parameter. For 
large values of l^l, the LEP 190 curve roughly follows the LEP 175 curve as discussed in the text. 
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